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TECHNICAL  REPORT  SUMMARY 


/  -^yThis  project  examined  the  contribution  of  scattering  to  the  attenu¬ 
ation  of  short  pulses  within  the  crust.  Coda  decay  and  excitation  for 
local  events  were  examined  at  Mammoth  Lakes  and  Morgan  Hill,  CallfdmUC^ 
Monticello,  S0utl/^jfeal4<wt^and^ New  Brunswick,  Canada,  in  the  frequency 

range  3-50  Hz.  For  short  times  (less  than  10  seconds),  the  total  tur- 

■  /  - 

bidlty  determined  from  coda  decay  was  about  0.1/km  for  all  regions, 

..  ■  s<>i>  j,  >  / 

applying  a  magnitude  bias  of  0.2  in^  if  10  km  of  such  material  is 
traversed.  Since  the  total  turbidity  Is  Independent  of  frequency. 
Implying  geometrical  scattering,  this  would  not  be  detectable  by  spec¬ 
tral  ratio  methods.  The  back scattering  turbidity  determined  from  coda 
excitation  at  short  times  indicates  strong  scattering  in  the  upper 
crust,  especially  for  frequencies  in  the  3-10  Hz  range.  At  times  longer 
than  10-15  seconds  for  the  codas  from  the  eastern  North  American 


regions,  Monticello  and  New  Brunswick,  the  coda  energy  appeared  to  be 
channeled  Into  a  horizontally  propagating  mode  such  as  Lg.  The  total 
turbidity  for  this  portion  of  the  coda  was  lower  than  for  the  short 
codas,  about  O.O^knT^lndl eating  less  scattering,  a  result  born  out  by 
the  back scattering  turbidity.  Codas  from  California,  however,  did  not 
show  this  phenomenon.  Indicating  either  that  this  mode  is  not  present  or 
that  It  Is  more  strongly  scattered.  This  result  indicates  that  attenua¬ 
tion  for  Lg  can  be  estimated  from  the  coda  after  10  seconds,  as  proposed 
by  other  workers. a  All  of  the  scattering  observed  In  this  study  occurs 
In  the  crust;  slnce^ln  situations  where  single  scattering  occurs  the 


back scattering  turbidity  is  about  1$  of  the  total  turbidity,  velocity 
fluctuations  may  be  responsible.  Multiple  scattering  occurred  In  the 
short  codas  at  frequencies  around  3  Hz. 


INFLUENCE  OF  SCATTERING  ON  Q  IN  THE  LITHOSPHERE 

By 

Anton  M.  Dainty,  Robert  M.  Duckworth,  and  An  Tie 

Introduction 

If  accurate  assessments  are  to  be  made  of  the  size  and  nature  of  a 
seismic  event  (explosion  or  earthquake),  we  must  have  a  good  under¬ 
standing  of  the  changes  in  amplitude  and  waveform  of  seismic  waves  as 
they  travel  from  the  source  to  a  more  or  less  distant  receiver.  Conven¬ 
tionally,  this  problem  has  been  divided  into  two  parts:  "geometrical 
spreading"  and  "attenuation".  These  terms  derive  from  a  combination  of 
theoretical  considerations  of  idealised  problems  (ray  theory  in  homo¬ 
geneous,  isotropic— at  least  piecewise— media)  and  some  simple  experi¬ 
mental  observations  showing  that  in  real  media  a  certain  amount  of 
seismic  energy  is  apparently  lost,  perhaps  to  heat,  compared  to  ray 
theory  calculations.  Anything  that  can  be  fitted  into  the  theory  of 
wave  propagation  In  Ideal,  elastic  media  Is  usually  called  "geometrical 
spreading",  while  anything  that  can  be  described  as  loss  of  energy  from 
the  seismic  wavefleld  Is  usually  called  "attenuation".  In  the  simplest 
cases,  the  effect  of  geometrical  spreading  on  amplitude  often  has  the 
form  of  a  frequency  Independent  power  law  (e.g.,  I/R  for  spherical 
spreading),  whereas  the  effect  of  attenuation  on  amplitude  is  described 
by  a  frequency  dependent  exponential  law  (exp[-irfR/QV],  where  Q  Is  the 
quality  factor  and  V  is  the  wave  velocity). 

The  real  media  encountered  In  the  earth  differ  from  Ideal  homo¬ 
geneous  elastic  media  In  several  ways.  One  Is  the  phenomenon  of  conver- 


slon  of  the  elastic  energy  of  the  seismic  wavefleld  to  heat  mentioned 
above.  Another  is  the  presence  of  Inhomogeneity  of  the  elastic  para¬ 
meters  and  the  density,  probably  on  all  linear  scales.  Some  of  the 
effects  of  the  large  scale  inhomogeneity  can  be  Included  in  geometrical 
spreading  explicitly— for  example,  the  depth  dependent  Inhomogeneity 
commonly  observed  In  the  earth.  However,  much  of  the  inhomogeneity  is 
too  complex  to  be  handled  by  deterministic  methods  and  has  to  be  con¬ 
sidered  as  producing  random  scattering,  whose  effects  can  only  be  calcu¬ 
lated  as  a  statistical  average.  In  considering  these  effects,  a  diffi¬ 
culty  emerges.  Whilst  conceptually  the  effect  of  scattering  is  geo¬ 
metrical  spreading,  since  energy  does  not  leave  the  total  seismic  wave- 
field,  operationally  the  effect  of  scattering  may  in  part  be  like  that 
of  attenuation.  For  example,  If  a  pulse  confined  in  time  propagates 
through  a  scattering  medium,  seismic  energy  will  be  scattered  out  of  the 
time  frame  containing  the  pulse.  If  amplitude  measurements  are  made 
only  In  the  time  frame  containing  the  pulse,  a  common  procedure,  the 
amplitude  will  appear  to  have  attenuated  with  an  exponential  law  as 
exp[-GR/2],  where  G  is  a  parameter  known  as  the  total  turbidity,  or  the 
total  scattering  cross-section  per  unit  volume. 

In  the  work  reported  on  here,  we  have  endeavored  to  find  the  atten¬ 
uation  effect  due  to  scattering  in  the  lithosphere  of  the  earth  at  high 
frequencies  (above  1  Hz)  In  the  sense  described  above.  Since  sources  of 
particular  Interest  (explosions)  and  receivers  are  located  near  the 
surface  of  the  earth,  all  seismic  waves  encounter  this  environment;  for 
regional  phases  this  Is  the  only  environment  encountered.  Scattering 
has  been  assessed  by  examination  of  the  coda  of  local  earthquakes. 
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Since  the  coda  consists  of  backscattered  energy,  we  have  separately 
assessed  the  backscatterlng  turbidity  (or  back scattering  cross-section 
per  unit  volume)  determined  from  the  size  of  the  coda  relative  to  the 
direct  S  wave  and  the  total  quality  factor  Q,  expressed  as  an  apparent 
total  turbidity,  determined  from  the  coda  decay  rate.  We  have  done  this 
for  two  regions  of  California— Mammoth  Lakes  and  Morgan  Hill— and  two 
regions  in  eastern  North  America— Monti  cello.  South  Carolina,  and  New 
Brunswick,  Canada. 

Theoretical  Sumnary 

In  this  section  we  will  discuss,  in  order,  the  characterisation  of 
scatterers  and  scattering,  attenuation  of  a  short  pulse  due  to  scat¬ 
tering,  and  the  theory  of  local  coda.  Part  of  this  material  Is  covered 
in  earlier  reports  submitted  under  this  grant  and  included  as  Appendices 
1  and  2— references  will  be  made  to  the  Appendices  as  appropriate. 

The  simplest  case  of  scattering  occurs  when  there  is  only  one 
scatterer  of  prescribed  shape.  Whilst  the  case  of  scattering  in  an 
elastic  medium  has  recently  been  treated  by  Wu  and  Aki  (1984),  many 
useful  ideas  about  scattering  can  be  Illustrated  by  simpler  acoustic 
cases.  Let  us  consider  a  sound  wave  in  a  uniform  acoustic  medium  of 
amplitude  Aj  and  frequency  f  Incident  on  a  body  of  different  acoustic 
properties  located  at  the  origin  of  spherical  polar  coordinate  R,  0,  <f>. 
At  large  distances  from  the  body  the  amplitude  of  the  scattered  sound 
wave  Ag  will  be  given  by  (Clay  and  Medwln,  1977): 


A?  -  A?  •  o(e,*,f)  •  (1/R2) 


(1) 
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This  indicates  that  far  away  from  the  scattering  body  the  body  appears 
as  a  point  source  of  magnitude  Aj  v,a(e,4>,f),  with  an  angular  dependence 
of  radiation  controlled  by  the  quantity  o(e,  <f>,f)  known  as  the 
differential  cross-section.  An  important  special  case  is  e  *  v 
(back scattering  in  the  strict  sense),  then  a  *  ab(f)»  the  backscattering 
cross-section.  Figure  1  shows  ab  for  a  rigid,  fixed  sphere  in  an 
acoustic  medium.  An  important  feature  of  Figure  1  is  the  division  of 
scattering  into  two  types.  At  low  frequencies,  the  backscattering  cross 
section  increases  strongly  (as  f*)  with  frequency  (Rayleigh  scattering) 
as  long  as  the  radius  of  the  sphere  is  appreciably  less  than  the  wave¬ 
length.  Geometrical  scattering  occurs  at  higher  frequencies  when  the 
wavelength  is  shorter  than  the  radius  of  the  sphere.  Under  these  condi¬ 
tions,  the  backscattering  cross-section  for  geometrical  scatter  is 
larger  than  for  Rayleigh  scatter,  but  independent  of  frequency.  Also, 
for  both  cases,  the  backscattering  cross-section  is  proportional  to  the 
physical  cross-section.  These  results  taken  together  indicate  that  if  a 
variety  of  scatterers  of  different  sizes  are  Insonified  by  a  sound  wave 
of  a  given  wavelength,  the  larger  scatterers  that  are  in  the  geometrical 
scattering  regime  will  dominate,  other  things  being  equal. 

If  there  are  N  scatterers  per  unit  volume,  the  differential  cross- 
section  per  unit  volume,  or  differential  turbidity,  is 

g(e,$,f)  ■  No( e, $,f )  (2) 

The  energy  scattered  in  direction  0,41  from  a  small  volume  aV  is  the 


Aj  *  Aj  g(M,f)  ‘  AV/R2 


(3) 


This  formulation  leads  to  a  second  method  of  describing  scatterers,  in 
terms  of  a  random  medium.  Chernov  (I960),  using  results  from  earlier 
workers,  showed  that  if  an  acoustic  medium  contained  velocity  or  density 
fluctuations  ("scatterers"),  only  a  knowledge  of  the  autocorrelation  of 
these  fluctuations  is  necessary  to  calculate  g(6,4>,f). 

To  apply  (3)  to  the  attenuation  of  a  short  pulse,  note  that 

scattered  energy  may  be  effectively  lost  from  the  pulse.  In  the 

discussion  that  follows,  the  pulse  will  be  discussed  in  terms  of  its 

Fourier  components.  Energy  deflected  at  small  angles  to  the  original 

direction  of  travel  (forward  scatter)  will  not  be  lost,  however;  Sato 

(1982)  suggests  that  most  energy  scattered  at  angles  of  30°  or  less  is 

included  in  the  original  pulse,  causing  fluctuations  of  amplitude  and 

phase.  The  total  loss  of  intensity  from  a  wave  encountering  a  volume 

AV  of  scatterers  can  be  found  by  enclosing  the  volume  in  a  sphere  of 

radius  R  and  integrating  the  scattered  intensity  given  by  (3)  over  the 

surface  of  the  sphere,  excluding  energy  scattered  at  angle  0  less  than 

e„  ~  30°,  as  discussed  above.  Here  0  ■  0  is  taken  as  the  direction  of 
c 

propagation  of  the  incident  wave.  Then  the  intensity  loss 


Al j  s  a(Aj )  =  -A i 


d4>  J  de  •  sins  *  g(0,<f>,f) 

e 

c 


=  -Ij  •  AV  .  G(f) 


(4) 


In  (4)  the  I  * s  represent  intensity  and  G(f)  is  known  as  the  total 
turbidity.  Using  (2),  we  may  also  define 


the  total  scattering  cross-section. 


If  AV  In  (4)  Is  taken  as  a  slab  consisting  of  a  unit  area  face  in 
the  plane  of  the  wavefront  and  a  side  Ax  perpendicular  to  the  wavefront 

Al{  -  -6(f)  •  Ij  •  Ax  (6) 

If  we  consider  the  wave  to  propagate  through  a  scattering  medium  that 
may  be  split  up  Into  such  slabs,  (6)  can  be  integrated  (Wu,  1982)  to 
obtain 


Ij  *  IQ  exp[-G(f)x]  (7) 

This  may  be  expressed  formally  In  terms  of  a  scattering  Q,  Qs,  If  Q$  is 
defined  by 

Ij  *  IQ  expC-2,fx/(QsC)]  (8) 

where  C  Is  the  sound  velocity.  Then 

1/QS  -  G(f)C/(2irf)  (9) 

If  the  medium  also  attenuates  anelastically,  parameterlsed  by  an 
intrinsic  Q,  Qj ,  then  the  total  Q  will  be  given  by  (Dainty,  1981) 

1/Q  -  1/Qj  +  1/QS  (10) 

-  1/Qj  +  G(f  )C/(2irf )  (11) 

The  frequency  dependence  expected  for  Q$  may  be  deduced  from  the 
results  presented  earlier  for  a  rigid.  Immovable  sphere.  For  low 
frequencies  Rayleigh  scattering  occurs;  generally,  G(f)  has  a  similar 
frequency  dependence  as  g(ir,f),  that  Is  to  say,  G(f)  Increases  with 
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frequency  as  f  leading  to  Q$  a  1/f  .  For  high  frequencies,  the 
wavelength  will  be  short  and  geometrical  scattering  may  occur.  Then 
G(f)  will  be  Independent  of  frequency  and  Q$  a  f.  Dainty  (1981)  has 
examined  two  cases  of  shear  wave  Q  In  which  this  may  be  true  for 
frequencies  above  1  Hz  In  the  lithosphere.  Dainty  (1984)  and  Vlu  (1982) 
discuss  the  frequency  dependence  expected  for  Q$  when  the  scattering 
medium  is  characterised  as  a  random  medium  rather  than  by  discrete 
scatterers.  Dainty  (1984;  Included  as  Appendix  2  of  this  report) 
demonstrated  that  if  there  Is  a  length  scale  associated  with  the  random 
medium,  and  discontinuities  within  the  medium,  then  at  high  frequencies 
(wavelength  shorter  than  the  length  scale)  the  frequency  dependence  of 
g(9,<|>,f)  and  G(f)  will  be  the  same  as  that  for  geometrical  scattering, 
l.e. ,  independent  of  frequency. 

The  theory  presented  so  far  Is  for  the  acoustic  case.  Wu  and  Aki 
(1984)  have  presented  results  for  the  elastic  case.  They  find  that.  In 
addition  to  P  to  P  and  S  to  S  scattering,  which  are  analogous  to 
acoustic  scattering,  there  Is  also  P  to  S  and  S  to  P  scattering  ("mode 
conversion").  However,  there  are  no  backscattered  mode  conversions 
(l.e.,  scattered  back  directly  along  the  ray  path).  P  to  P  and  S  to  S 
scattering  Is  similar  to  acoustic  scattering  If  the  random  medium 
consists  of  velocity  fluctuations  In  that,  at  high  frequencies,  most  of 
the  energy  Is  scattered  forward.  However,  if  the  random  medium  consists 
of  Impedance  fluctuations,  most  of  the  energy  Is  backscattered  at  high 
frequencies. 

In  this  report,  we  shall  be  dealing  with  scattered  seismic  waves 
that  form  the  coda  of  local  earthquakes.  These  waves  arrive  after  the 


the  total  scattering  cross-section. 

If  aV  in  (4)  is  taken  as  a  slab  consisting  of  a  unit  area  face  In 
the  plane  of  the  wavefront  and  a  side  Ax  perpendicular  to  the  wavefront 

Alj  -  -G(f)  •  Ij  •  Ax  (6 

If  we  consider  the  wave  to  propagate  through  a  scattering  medium  that 
may  be  split  up  Into  such  slabs,  (6)  can  be  integrated  (Wu,  1982)  to 
obtain 

Ij  »  IQ  exp[-G(f)x]  (7 

This  may  be  expressed  formally  In  terms  of  a  scattering  Q,  Q$,  If  Qs  Is 
defined  by 

Ij  •  IQ  exp[-2wfx/(QsC)]  (8 

where  C  is  the  sound  velocity.  Then 

1/QS  -  G(f  }C/(2irf )  (9 

If  the  medium  also  attenuates  anelastically,  parameterised  by  an 
intrinsic  Q,  Qj,  then  the  total  Q  will  be  given  by  (Dainty,  1981) 

1/Q  -  1/Qj  +  1/QS  (10 

-  1/Qj  +  G(f )C/(2wf )  (11 

The  frequency  dependence  expected  for  Q$  may  be  deduced  from  the 
results  presented  earlier  for  a  rigid.  Immovable  sphere.  For  low 
frequencies  Rayleigh  scattering  occurs;  generally,  G(f)  has  a  similar 
frequency  dependence  as  g(ir,f),  that  is  to  say,  G(f)  increases  with 


4  3 

frequency  as  f  leading  to  Qs  a  1/f  .  For  high  frequencies,  the 
wavelength  will  be  short  and  geometrical  scattering  may  occur.  Then 
G(f)  will  be  Independent  of  frequency  and  Qg  o  f.  Dainty  (1981)  has 
examined  two  cases  of  shear  wave  Q  in  which  this  may  be  true  for 
frequencies  above  1  Hz  In  the  lithosphere.  Dainty  (1984)  and  Wu  (1982) 
discuss  the  frequency  dependence  expected  for  when  the  scattering 
medium  Is  characterised  as  a  random  medium  rather  than  by  discrete 
scatterers.  Dainty  (1984;  Included  as  Appendix  2  of  this  report) 
demonstrated  that  If  there  Is  a  length  scale  associated  with  the  random 
medium,  and  discontinuities  within  the  medium,  then  at  high  frequencies 
(wavelength  shorter  than  the  length  scale)  the  frequency  dependence  of 
g(M,f)  and  G(f)  will  be  the  same  as  that  for  geometrical  scattering, 
l.e. ,  Independent  of  frequency. 

The  theory  presented  so  far  Is  for  the  acoustic  case.  Wu  and  Aki 
(1984)  have  presented  results  for  the  elastic  case.  They  find  that.  In 
addition  to  P  to  P  and  S  to  S  scattering,  which  are  analogous  to 
acoustic  scattering,  there  Is  also  P  to  S  and  S  to  P  scattering  ("mode 
conversion").  However,  there  are  no  backscattered  mode  conversions 
(l.e.,  scattered  back  directly  along  the  ray  path).  P  to  P  and  S  to  S 
scattering  Is  similar  to  acoustic  scattering  if  the  random  medium 
consists  of  velocity  fluctuations  In  that,  at  high  frequencies,  most  of 
the  energy  Is  scattered  forward.  However,  if  the  random  medium  consists 
of  Impedance  fluctuations,  most  of  the  energy  is  backscattered  at  high 
frequencies. 


In  this  report,  we  shall  be  dealing  with  scattered  seismic  waves 
that  form  the  coda  of  local  earthquakes.  These  waves  arrive  after  the 


latest  direct  arrivals,  are  Incoherent  between  components,  and  are 
believed  to  be  either  backscattered  surface  waves  or  backscattered  S 
waves  (Akl  and  Chouet,  1975).  The  geometry  of  the  backscattered  paths 
considered  in  this  model  is  shown  In  Figure  2.  If  times  are  long  enough 
the  source  and  receiver  may  be  considered  to  be  coincident.  The  source 
pulse  is  considered  to  be  short  and  equivalent  to  the  S  wave  pulse. 

Only  single  scattering  Is  considered.  The  theory  for  this  case  Is  given 
In  Akl  and  Chouet  (1975)  and  Appendix  1  attached  to  this  report.  For  S 
to  S  scattering  In  a  medium  of  average  velocity  V  (three  dimensional 
spreading),  the  power  spectrum  as  a  function  of  time  t  Is  given  by: 

P(f »t)  »  Ab(f)t“2  exp[-2»ft/Q]  (12) 

Ab(f)  -  S2(f)  •  exp[2wfts/Q]  •  8*V  •  g  (w,f)  •  t2  (13) 

S(f )  Is  the  Fourier  transform  of  the  S  wave  pulse,  t$  Is  the  travel  time 
of  the  direct  S  wave,  g(*,f)  Is  the  S  to  S  back scattering  turbidity  and 
Q  is  given  by  (10)  or  (11).  For  the  case  of  surface  wave  scattering 
(two  dimensional  spreading), 

P(f.t)  -  A$(f )t-1  exp[-2*ft/Q']  (14) 

A$(f)  -  CS'(f)]2  •  exp[2*ftj/Q']  •  ZwV  .  g‘(„,f)  .  t;  (15) 

The  primed  variables  In  (14)  and  (15)  have  the  same  definitions  as 
before,  but  for  surface  waves  rather  than  for  S  waves.  The  results 
discussed  In  this  report  will  be  Interpreted  using  (12)  and  (14). 


The  basic  data  base  consists  of  digital  seismograms  of  local 
earthquakes.  Five  sets  of  seismograms  were  available.  Four  of  these 
sets  were  taken  by  the  United  States  Geological  Survey,  or  using  their 
equipment,  and  consist  of  sampled  time  series  of  three  component  data 
with  a  sampling  Interval  of  0.005  sec.  These  sets  were  taken  at  Mammoth 
Lakes,  California;  Morgan  Hill,  California;  Montlcello,  South  Carolina; 
and  New  Brunswick,  Canada.  The  fifth  set  consists  of  digital  seismo¬ 
grams  taken  by  the  Eastern  Canadian  Network,  McKendrlck  Lake  Station, 
vertical  component  only,  recorded  at  a  sampling  Interval  of  60  samples 
per  second.  Procedures  for  analysing  the  data  according  to  equations 
(12)  or  (14)  have  been  presented  In  Duckworth  and  Dainty  (1983),  In¬ 
cluded  In  this  report  as  Appendix  1.  Briefly,  one  component  of  the 
seismogram  Is  Fourier  transformed  In  a  sliding  window  and  the  average 
Fourier  amplitude  as  a  function  of  time  determined.  Equation  (12)  or 
(14)  is  then  fitted  to  the  coda  at  times  greater  than  twice  the  S  wave 
travel  time,  estimated,  together  with  the  origin  time,  from  the  S-P 
times.  The  S  wave  Fourier  amplitude  Is  also  determined  by  taking  the 
average  of  the  two  nearest  windows;  the  window  length  (usually  1.28  sec) 
was  chosen  to  correspond  roughly  to  the  S  pulse  length.  Then  the  two 
unknown  parameters,  Q  and  g(w,f)  (or  their  primed  equivalents),  were 
determined. 

Preliminary  results  obtained  by  Duckworth  and  Dainty  (1983) 
Indicated  that  Q  was  approximately  proportional  to  frequency  over  the 
frequency  range  3-50  Hz.  This  suggests  attenuation  primarily  by 
geometrical  scattering  (see  Dainty,  1984,  Appendix  2).  Accordingly,  to 


compare  Q  and  g(w,f),  a  quantity  called  the  apparent  total  turbidity  is 
calculated  from 

Ga(f)  *  2wf/(QV)  (16) 

V  in  this  report  is  taken  to  be  3.2  km/sec,  typical  of  shear/surface 
wave  velocities  In  the  upper  crust.  If  anelastic  attenuation  can  be 
neglected.  6  (f)  will  be  equal  to  G(f).  Otherwise, 

a 

Ga(f)  >  6(f)  >  g(w.f)  (17) 

The  last  Inequality  follows  from  (3)  and  (4). 

In  attempting  to  fit  codas  from  Monti  cello.  South  Carolina,  it  was 
found  that  a  single  value  of  A^  (or  A$)  and  G(»,f)  (or  G'(ir,f))  did  not 
provide  a  good  fit  to  codas  of  total  duration  longer  than  10  seconds; 
there  appeared  to  be  a  change  In  the  decay  rate  of  the  coda  at  this  time 
(see  Appendix  1  for  an  example).  For  these  codas,  separate  fits  were 
made  of  the  coda  before  10  seconds  and  the  coda  after  10  seconds.  A 
similar  effect  was  noted  on  the  records  from  the  Eastern  Canadian 
Network  station  with  the  change  in  decay  rate  occurring  at  about  15 
seconds.  Since  the  McKendrick  Lake  station  was  25  km  from  the  after¬ 
shock  area,  meaning  that  the  direct  S  travel  time  was  about  7.5  seconds. 
It  was  not  possible  to  fit  the  early  coda  and  maintain  the  rule  that 
only  coda  at  times  greater  than  twice  the  direct  S  travel  time  Is  used 
for  the  fit.  However,  the  results  from  McKendrick  Lake  may  be  compared 
to  results  from  the  United  States  Geological  Survey  New  Brunswick  seis¬ 
mograms,  which  are  In  the  same  area  and  from  the  same  aftershock 


sequence,  but  much  closer.  The  U.S.G.S.  seismograms  are  too  short  to 
cover  the  time  frame  beyond  15  seconds. 


In  fitting  (12)  or  (14),  only  one  component  of  data  has  been  used, 
usually  the  vertical.  Properly,  the  sum  of  the  squared  amplitude 
spectra  for  three  components  In  each  window  should  be  used.  However, 
three  component  data  was  not  always  available,  so  single  component  data 
was  analysed  so  that  results  could  be  compared.  This  should  not  affect 
the  measurement  of  G  (f),  since  the  decay  rate  should  be  the  same  for 

d 

all  components  provided  that  one  Is  not  examining  different  types  of 
waves  on  different  components,  but  the  measurement  of  g(ir,f)  could  be 
affected.  To  test  this,  and  to  examine  the  effect  of  different 
components  on  estimation  of  G_(f),  analyses  of  three  component 

d 

measurements  carried  out  separately  for  each  component  are  presented  in 
Figure  3  for  a  Morgan  Hill  event  and  in  Figure  4  for  a  New  Brunswick 
event,  using  U.S.G.S.  data  In  the  frequency  range  3-50  Hz.  In  both 
cases,  the  results  for  G  (f)  are  quite  consistent  for  the  three 
components,  with  the  exception  of  the  vertical  in  New  Brunswick  for  the 
frequencies  3  and  6  Hz.  The  results  for  g(*,f)  are  also.  In  the  main, 
quite  consistent.  This  Indicates  that  sufficiently  accurate  results  for 
our  purposes  may  be  found  by  analysing  only  a  single  component  of  the 
data;  this  considerably  Increases  the  amount  of  data  that  may  be 
considered. 

An  Implicit  assumption  made  in  this  analysis  Is  that  the  parameters 
G  (f)  and  g(»,f)  derived  are  characteristic  of  the  region  In  which  the 
measurements  are  taken  and  not  the  station  or  the  event.  This 
assumption  was  tested  at  the  Mammoth  Lakes,  California,  site,  and  the 


lew  Brunswick,  Canada.  Symbols  as  for  Figure 


results  are  shown  in  Figures  5  and  6.  While  there  is  considerable 
variation,  the  values  are  consistent,  within  an  order  of 
magnitude. 

Typical  results  for  Morgan  Hill,  California,  and  New  Brunswick, 
Canada  (codas  less  than  15  seconds  long),  have  already  been  shown 
(Figures  3  and  4),  using  equation  (14)  (cylindrical  spreading).  Figure 
7  shows  typical  results  from  New  Brunswick,  Canada  (codas  more  than  15 
seconds  long),  again  using  equation  (14).  Figures  8,  9,  and  10,  taken 
from  Dainty  (1984),  show  averaged  results  for,  respectively,  Monticello, 
South  Carolina  (codas  less  than  10  seconds  long);  Monticello,  South 
Carolina  (codas  more  than  10  seconds  long);  and  Mammoth  Lakes, 
California. 

Discussion 

From  Figures  3,  4,  7,  8,  9,  and  10,  the  codas  examined  here  may  be 
divided  Into  two  cases.  Figures  7  and  9  show  the  results  from  long 
codas  in  eastern  North  America  (Monticello,  South  Carolina,  and  New 
Brunswick,  Canada),  a  non-tectonlc  area.  The  values  of  Gff)  are  mostly 

a 

between  0.005  and  0.05  km”*,  and  the  values  of  g(*,f)  are  between  0.0001 
and  0.002  km”1,  substantially  less  than  most  of  the  corresponding  values 
In  Figures  3,  4,  8,  and  10.  Also,  in  fitting  (12)  and  (14)  to  the  data, 
negative  values  of  Q  (and  G  (f ) )  were  often  obtained  for  (12)  (spherical 
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spreading).  This  Is  Indicated  In  Figure  7  (Monticello  data)  for 
frequencies  of  6,  9,  and  12  Hz;  for  the  New  Brunswick  data,  Q  was 
negative  for  3  Hz,  Infinite  for  6  Hz,  and  very  high  for  the  two  higher 
frequencies.  From  (12),  a  negative  Q  means  that  the  coda  power  spectrum 
decays  less  rapidly  than  t  .  This  probably  means  that  these  long 


closed  squares,  TOM.  Results  above  30  Hz  probably  not  reliable. 


spreading  assumed,  closed  symbols,  cylindrical  spreading  assumed.  Circles 


than  10  seconds  long.  Symbols  as  for  Figure  8  . 
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codas  consist  of  energy  that  is  at  least  partially  "channeled" 
horizontally,  like  an  Lg  wave.  This  in  turn  implies,  that  measurements 
on  such  codas  (long  means  longer  than  10  or  15  seconds;  this  is  in  fact 
the  usual  situation)  would  probably  give  a  good  estimate  of  Lg  Q 
(Mitchell  and  Nuttll,  1983;  Pomeroy  and  Peseckis,  1984).  Also,  it 
confirms  the  importance  of  scattering  In  Lg  propagation  (Gupta  and 
Blandford,  1983). 

The  short  codas  (less  than  10  or  15  seconds  long)  in  eastern  North 
America  (New  Brunswick,  Figure  4;  Montlcello,  Figure  8)  and  all  the 
codas,  long  or  short,  in  California  (Morgan  Hill,  Figure  3;  Mammoth 
Lakes,  Figure  10)  have  higher  rates  of  decay,  as  shown  by  values  of 
G  (f)  of  around  0.1  km"1  and  generally  higher  values  of  g(ir,f),  between 

a 

0.005  and  0.1,  than  the  long  codas. 

The  separation  of  long  and  short  codas  may  be  explained  in  two 
ways.  One  possibility  is  that,  in  the  case  of  the  eastern  North  America 
codas,  there  are  two  modes  of  propagation,  one  of  them  ("short  codas") 
more  highly  scattered  (higher  G_ (f )  and  g(w,f))  than  the  other.  The 

Q  ““ 

more  highly  scattered  mode  will  dominate  at  short  times,  since  g(ir,f), 
which  controls  the  excitation  of  the  coda,  is  larger,  but  the  faster 
decay  of  this  mode  means  that  the  mode  with  the  smaller  value  of  G  (f ) 

Q 

(higher  Q,  by  (16))  will  dominate  at  long  times.  In  California,  either 
the  second  mode  Is  not  present  or  its  scattering  characteristics  are 
similar  to  those  of  the  first  mode.  If  the  Identification  of  the  long 
coda  mode  In  eastern  North  America  with  waves  similar  to  Lg  is  correct, 
the  absence  or  greater  scattering  of  such  waves  in  California  Implies 


that  Lg  would  not  propagate  as  well  In  California  as  In  eastern  North 
America,  In  agreement  with  observation  (Nuttll,  1973). 

A  second  possibility  Is  that  the  change  In  coda  characteristics  at 
Montlcello  and  New  Brunswick  Is  due  to  a  change  from  single  scattering 
to  multiple  scattering  at  10-15  seconds  time.  Sato  (1977)  gives 

t  >  1/(GV)  (18) 

as  a  necessary  but  not  a  sufficient  condition  for  multiple  scattering. 

If  6  -v  0.1  km”1,  a  typical  value  for  the  short  codas,  then  (18)  is  met 
for  t  >  3  sec.  Also,  for  the  extreme  case  of  very  strong  scattering, 
seismic  energy  In  the  coda  will  diffuse  (Dainty  and  Toksoz,  1977),  and 
an  appropriate  formula  for  the  coda  power  spectrum  in  a  three- 
dimensional  medium  Is 

P(f,t)  -  E(f)t’3/2  exp[-3R2G/(4Vt)]  exp[-2»ft/Q1 ]  (19) 

-  E(f)t’3^2  exp[-2irft/Qj  ]  ,  R2/t  -►  0  (20) 

E(f)  *  S2(f)  ;  / 27G3  /  (8  / w3V3  )  (21) 

R  Is  the  source-receiver  distance,  taken  to  be  small  In  (20).  Like  (12) 

and  (14),  (20)  has  the  form  of  an  Inverse  power  law  times  an 

exponential,  but  the  exponential  contains  only  the  Intrinsic  attenuation 

Q.|,  not  Q  defined  by  (11).  Accordingly,  a  change  In  decay  rate  to  lower 

values  would  be  expected  If  the  coda  changed  from  a  single  scattering 

phenomenon  to  a  multiple  scattering  phenomenon.  Finally,  we  note  that 

for  short  codas  at  frequencies  less  than  10  Hz,  relation  (17)  is  often 

not  met,  at  the  same  time  that  G  (f)  Is  decreasing  relative  to  Its  value 
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at  higher  frequencies.  This  Is  a  strong  Indication  that  multiple 
scattering  Is  occurring  In  the  short  codas,  and  that  the  values  of  G  (f) 

u 

obtained  by  fitting  (12)  or  (14)  are  too  low. 

In  spite  of  these  Indications  that  multiple  scattering  Is  occurring 
In  the  short  codas,  the  suggestion  that  the  change  from  short  to  long 
codas  reflects  a  change  from  single  to  multiple  scattering  Is  rejected 
for  the  following  reasons.  First,  all  of  the  arguments  In  favor  of 
multiple  scattering  In  the  eastern  United  States  apply  with  equal  force 
to  California,  yet  the  change  In  coda  Is  not  observed  there.  Thus, 
there  must  be  regional  differences.  Furthermore,  If  condition  (18)  is 
applied  with  G  ^  0.01  km'*,  a  typical  value  for  the  Monticello  codas, 
then  single  scattering  would  be  expected  for  t  <  30  sec,  which  is  the 
typical  length  of  one  of  these  codas.  For  these  Monticello  codas, 
condition  (17)  Is  also  met,  although  not  at  3  Hz  for  the  New  Brunswick 
codas. 

The  nature  of  the  scatterers  and  the  contribution  of  scattering  to 
attenuation  may  be  partially  answered  from  the  data  presented  here.  For 
the  long  codas  at  Monticello,  there  Is  a  good  correlation  between  G  (f) 
and  g(w,f),  particularly  around  10  Hz.  For  all  codas  except  the  long 
Monticello  codas,  there  Is  a  good  evidence  that  scattering  is  strong  for 
frequencies  below  about  10  Hz,  which  would  Indicate  a  substantial 
contribution  to  attenuation  through  (10)  or  (11).  In  cases  below  10  Hz 
where  multiple  scattering  Is  believed  to  be  occurring,  G  (f)  declines 
relative  to  higher  frequencies,  again  Indicating  that  a  substantial  part 
of  the  attenuation  Is  due  to  scattering,  from  (11),  (16),  and  (20). 
Finally,  as  pointed  out  earlier,  G  (f)  Is  approximately  constant  with 


frequency.  In  the  main,  implying  Q  proportional  to  frequency;  this  is. 
characteristic  of  attenuation  by  geometrical  scattering.  For  these 
reasons,  we  estimate  that  at  least  50%  of  the  attenuation  observed  in 
the  crust  Is  due  to  scattering. 

The  nature  and  location  of  the  scatterers  may  be  partly  determined 
from  the  data  presented  here.  In  cases  where  single  scattering  Is 
believed  to  be  occurring,  l.e.,  long  Monti  cello  codas,  and  other  codas 
at  high  frequencies  ( >10  Hz),  g(*,f)  Is  about  1%  of  G  (f);  this  is 

d 

typical  of  velocity  fluctuations  In  a  three  dimensional  medium  (Uu  and 
Aki,  1984;  Dainty,  1984,  Appendix  2),  although  It  is  not  clear  what 
would  happen  in  a  two  dimensional  case.  Since  the  events  range  from  0 
to  10  km  depth,  the  short  codas  must  be  at  least  partially  due  to 
scatterers  In  the  upper  crust,  since  these  codas  start  Immediately  after 
the  direct  S  pulse.  The  long  codas  start  at  10-15  seconds,  which 
implies  a  mid  to  upper  crust  propagation,  since  the  scattered  wave  has 
to  get  to  the  scatterer  from  the  source  and  back  in  10-15  seconds  at 
^3.2  km/sec. 

The  principal  remaining  question  with  the  codas  studied  here  is 
whether  horizontal  channeling  of  energy  (two  dimensional  spreading, 

(14))  as  opposed  to  spherical  spreading  (12)  is  the  mode  of  energy 
propagation.  For  the  long  codas,  evidence  has  been  offered  that 
horizontal  channeling  Is  occurring;  for  the  short  codas,  the  question  is 
open.  Note  that  for  horizontal  channeling,  more  than  surface  waves, 
sensu  strlctu  must  be  Involved,  since  the  skin  depth  of  surface  waves 
at,  say,  10  Hz  Is  a  few  hundred  meters,  less  than  the  focal  depth  of 


many  of  the  events.  Generally,  the  decomposition  of  codas  into  waves  of 
different  types  Is  a  major  remaining  problem. 

What  is  the  significance  of  these  results  for  the  estimation  of 
seismic  yield?  First,  they  illustrate  that  scattering  occurs  in  the 
upper  part  of  the  crust,  where  the  source  and  the  receiver  are  located, 
and  this  scattering  has  a  significant  effect  on  attenuation.  For 
example.  If  G  (f)  is  0.1,  then  on  travelling  through  10  km  of  such 

Q 

material,  the  energy  in  a  short  pulse  will  be  down  by  1/e.  This  is 
equivalent  to  about  0.2  on  the  m^  magnitude  scale,  which  uses 
measurements  at  frequencies  (M.  Hz)  close  to  those  investigated  here  (3- 
50  Hz).  If  the  scattering  effect  varied  from  place  to  place,  this  could 
introduce  bias  of  the  order  of  tenths  of  a  magnitude  point,  which  is 
significant.  Further,  if  G  (f)  Is  independent  of  frequency  over  some 

a 

frequency  range,  as  has  been  observed  in  this  study,  then  the 
attenuation  effect  could  not  be  detected  by  spectral  ratio  methods 
within  that  frequency  range. 

This  study  also  Indicates  that  the  attenuation  of  the  regional 
phase  Lg  should  vary  from  region  to  region,  being  low  In  regions  wimilar 
to  eastern  North  America,  a  non-tectonlc  region,  and  high  In  regions 
similar  to  California,  a  tectonic  region,  in  accord  with  observations 
(Nuttli,  1973).  A  good  estimate  of  this  attenuation  should  be 
obtainable  from  the  coda  after  10  seconds  (Mitchell  and  Nuttli,  1983). 

Three  questions  for  further  Investigation  suggest  themselves:  Can 
the  scattering  characteristics  derived  for  a  region  from  local  codas  be 
found  from  teleseismic  observations,  such  as  teleseismlc  P  codas?  Could 
a  region  be  "calibrated"  in  terms  of  the  scattering  attenuation  of 


%  \  * 
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teleseismlc  phases  by  observations  of  local  coda?  What  other  effects 
does  scattering  have  on  the  waveforms  of  short  pulses? 
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OBSERVATIONS  OF  COOA  Q  FOR  THE  CRUST 
NEAR  MAMMOTH  LAKES,  CALIFORNIA,  AND  MONTICELLO,  SOUTH  CAROLINA 

By  Robert  M.  Duckworth  and  Anton  M.  Dainty 

Introduction 

The  seismic  Q  of  a  material  is  a  measure  of  the  degree  to  which 
it  attenuates  the  harmonic  energy  content  of  seismic  waves  passing 
through  it.  Seismic  energy  released  by  earthquakes  and  underground 
explosions  propagates  in  the  earth  primarily  in  the  form  of  elastic 
waves.  To  the  extent  that  earth  materials  are  not  perfectly  elastic 
the  elastic  waves  will  be  attenuated.  The  energy  at  an  initial  time, 
Eq,  and  the  energy  E  after  elapsed  time  t  are  related  to  the  Q  by 

E  =  Eq  exp(-2irft/Q)  (1) 

where  f  is  the  frequency  of  the  wave.  Attenuation  is  proportional  to 
the  distance  a  wave  travels,  x,  since  t  =  x/B  where  B  is  the  speed  at 
which  the  wave  propagates.  Therefore,  knowledge  of  Q  should  be  useful 
in  predicting  the  expected  amplitudes  of  seismic  waves  for  yield 
estimation  and  other  applications.  However,  because  earth  materials 
are  inhomogeneous  a  prediction  of  seismic  wave  amplitude  is 
complicated. 

Inhomogeneities  lead  to  reflections  and  refractions,  or 
scattering,  of  seismic  waves.  Scattering  by  these  inhomogeneities  is 
an  important  part  of  the  seismic  trace  recorded  as  a  function  of  time. 
The  first  part  of  a  seismogram  of  a  local  (within  a  few  hundred  km  of 


the  seismic  recording  station)  earthquake  is  the  Primary  or  P-wave. 
This  is  a  congressional  wave  and  travels  at  speeds  up  to  about  8  km/s 
for  local  earthquakes.  A  second  major  feature  is  the  Secondary  or 
S-wave  which  is  a  shear  wave  and  travels  at  speeds  of  about 
P-speed/^T.  Between  the  P  and  S  arrivals  and  following  the  S  arrival 
additional  waves  arrive.  Their  amplitude  decays  in  a  generally 
exponential  manner.  Aki  (1969)  proposed  that  these  waves  are  the 
result  of  scattering  of  surface  waves  by  surface  inhomogeneities.  He 
gave  them  the  name  Coda  [Latin:  tail]  waves.  Later*  Aki  and  Chouet 
(1975)  considered  body  wave  scattering.  Also,  scattering  affects  the 
attenuation  of  seismic  waves  in  two  ways.  Energy  is  scattered  out  of 
the  primary  wave  decreasing  the  apparent  Q.  Some  of  this  scattered 
energy  arrives  at  the  receiver  in  the  form  of  coda  waves.  The 
amplitude  of  these  depends  upon  the  scattering  mechanism, 
characterised  by  the  turbidity  as  described  below,  and  the  inelastic 
Q. 

This  report  presents  estimates  of  the  apparent  Q  from  the  coda  of 
local  earthquakes  in  two  different  geologic  regions.  Attempts  are 
made  to  separate  the  effects  of  scattering  and  inelasticity  in  order 
to  determine  the  turbidity  and  inelastic  Q  for  each  region.  The 
excitation  of  the  coda  has  also  been  considered. 

Theory 


Coda  waves  are  generally  considered  to  be  backscattered  energy. 
A  model  from  Aki  and  Chouet  (1975)  of  coda  waves  as  single,  S-wave  to 


S-wave  scatter  from  a  random  distribution  of  scatterers  is  used  in 
this  report. 

Observations  show  the  coda  to  be  independent  of  source-receiver 
distance,  R,  for  local  coda  waves  (Aki,  1969;  Aki  and  Chouet,  1975; 
Rautian  and  Khalturin,  1978).  Sato  (1977)  lends  support  to  this 
observation  with  a  model  for  single  isotropic  scatter.  He  shows  the 
time  dependence  of  the  mean  energy  density  of  the  coda  near  the  source 
to  be  nearly  independent  of  source-receiver  distance  R  for  R  <<  Bt/2. 
B  is  the  velocity  of  S-waves  in  the  medium  and  t  is  their  travel  time 
measured  from  the  origin  time.  Because  the  coda  shape  is  nearly 
independent  of  source  receiver  distance  for  local  coda  waves,  the 
model  is  simplified  to  the  case  of  coincident  source  and  receiver. 

Knopoff  and  Hudson  (1967)  modeled  scatter  from  small  perturba¬ 
tions  in  an  elastic  medium.  They  found  the  amplitude  for  S  to  S 
scatter  to  dominate  that  of  P  to  P  scatter  by  (a/B)1,  for  geometric 
(high  frequency)  and  (o/B)^  for  Rayleigh  (low  frequency)  scatter  and  P 
to  S  and  S  to  P  scatter  to  be  “vanishingly  small"  for  geometric 
scatter,  a  and  B  represent  the  speeds  of  P  and  S  waves  respectively. 
For  the  crust  o  *  /T  B  Is  a  good  approximation  and  the  amplitude  of  S 
to  S  scatter  Is  greater  than  P  to  P  scatter  by  a  factor  of  3  for 
geometric  scatter.  The  observed  character  of  coda  waves  also  supports 
S  to  S  dominance  (Aki,  1980). 

Before  presenting  a  model  of  coda  waves  the  effects  of  attenua¬ 
tion  will  be  considered.  Plane  waves  traveling  in  a  medium  with  a 
random  distribution  of  scatterers  will  be  attenuated.  Let  EQ  repre¬ 
sent  the  average  energy  of  the  waves  before  scattering  and  E  represent 
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the  energy  of  plane  waves  traveling  in  the  same  direction  as  the 
primary  waves  after  scattering.  We  assume  the  primary  cause  of 
attenuation  by  scattering  will  be  backscatter.  For  (EQ-E)  «  EQ  over 
a  distance  equal  to  the  average  separation  of  "backscatterers",  the 
attenuation  of  primary  waves  is  given  by 

E/EQ  =  exp (-xPob).  (2) 

The  backscatter  cross  section  represents  the  average  of  the  product 
of  the  backscatter  function  and  the  projected  area  of  the  scatterers, 
or  the  average  apparent  area  of  the  scatterers.  p  is  the  number  of 
scatterers  per  unit  volume.  As  a ^  and  p  are  unknown  for  the  crust, 
their  product  is  often  replaced  by  a  scattering  coefficient,  g,  called 
the  turbidity.  The  turbidity  represents  the  fractional  loss  of  energy 
due  to  scatter  per  unit  distance  traveled.  The  term  was  introduced  by 
Chernov  (1960)  and  first  used  in  seismic  literature  by  Galkin  et  aK 
(1970).  They  employ  it  to  describe  their  observed  high  frequency 
spatial  variation  of  seismic  amplitudes  as  a  function  of  distance  from 
the  source.  Equation  (2)  may  also  be  applied  to  the  attenuation  a 
plane  wave  suffers  when  passing  through  a  medium  where  scattering  is 
occurring  as  a  result  of  random  fluctuations  in  the  medium  properties 
rather  than  discrete  "backscatterers".  In  this  case,  g  is  a  function 
of  the  statistical  properties  of  the  mediun  and  the  wavelength  (Wu, 
1982). 

In  the  earth  medium  there  will  also  be  attenuation  resulting  from 
inelasticity.  This  is  attenuation  that  occurs  as  elastic  energy  Is 
converted  to  heat.  For  (E  -E)  «  E.  over  a  wavelength,  the  inelastic 


(3) 


E/Eq  -  exp{-2irfx/BQ). 

Q  is  the  quality  factor,  B  is  the  speed  of  shear  waves,  x  is  the 
distance  they  travel  in  the  medium  and  f  is  frequency.  The  inelastic 
attenuation  mechanism  is  not  well  understood.  It  is  generally 
observed  that  Q  is  independent  of  frequency  over  the  range  of  seismic 
frequencies  studied.  This  is  in  agreement  with  the  results  of  Knopoff 
(1964)  who  studied  Q  in  the  laboratory  and  found  it  independent  of 
frequency  for  dry  materials.  These  two  attenuation  effects  can  be 
combined.  In  terms  of  the  time  the  energy  spends  in  the  medium 
t  =  x/B  and  for  constant  B, 

E/EQ  *  exp(-2uft/Q)  .  exp(-Btg)  (4a) 

*  exp[-t(2irf/QT)].  (4b) 

Qt  represents  the  combined  effects  of  scattering  and  Inelastic 
attenuation. 


1/Qt  =  1/Q.  +  gB/2irf  ,  (5) 

where  Is  the  Inelastic  Q  (Dainty,  1981).  Dainty  and  Toksoz  (1981) 
suggest  that  for  a  single  scatter  theory  to  apply  2xf/Q^  »  Bg  is  a 
necessary  condition,  as  otherwise  there  would  be  sufficient  energy  for 
higher  order  scattering. 

The  power  spectrum  of  the  coda  as  a  function  of  time  under  the 
conditions  discussed  above  has  been  derived  by  Aki  and  Chouet  (1975). 
If  the  spectrum  of  the  source  is  ^(f ( rQ)  at  a  reference  distance  r0 


from  the  source,  the  coda  power  spectrum  is 

P(f|t)  «  st"2"1  exp(-2irft/QT)  (6) 

where 

s  »  |*(f|r0)|2  4mirr2m(B/2)1-2m  .  g  (7) 

s  represents  the  source  tern,  the  turbidity  g  *  pob,  and  m  «  0.5  for 
cylindrical  (surface  wave)  or  m  =  1.0  for  spherical  (body  wave) 
spreading.  The  source  term,  s,  will  in  general  depend  on  the 
scatterers  through  its  dependence  on  g.  The  time  t  is  measured  from 
the  origin  time  of  the  event. 

The  average  coda  amplitude  at  a  particular  frequency  and  time 
A(f|t)  is  related  to  the  power  spectrum. 

A(f|t)  *  [P(f jt)Af]1/2  (8) 

Af  Is  the  band  width  considered  about  f.  Qj  can  therefore  be  measured 
from  the  coda  wave  amplitude  attenuation  using  equation  (6).  Sato 
(1977)  points  out  the  single  scatter  model  is  only  applicable  for  a 
coda  duration  t  <  1/gB.  Longer  times  imply  a  significant  contribution 
of  higher  order  scattering  which  effects  the  return  of  scattered 
energy  to  the  later  parts  of  the  coda.  The  decomposition  of  QT  to 
find  Qj  and  g  depends  upon  assumptions  about  their  frequency  depen¬ 
dence.  This  and  the  applicability  of  the  single  scatter  model  will  be 
considered  in  the  section  on  results. 
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Data 

Two  sets  of  earthquake  data  are  available  for  analysis  courtesy 

of  Paul  Spudich  (U.S.  Geological  Survey  at  Menlo  Park,  California). 

Both  sets  contain  digital  records  of  earthquakes  and  their  locations. 

Program  Hypoinverse  by  Klein  (1978)  is  the  method  of  location. 

Locations  and  primary  S-wave  travel  times  of  the  earthquakes  selected 

for  analysis  have  been  read  from  the  U.S.G.S.  tapes  (Tables  1  and  2). 

The  first  set  is  from  a  series  of  events  associated  with  the  filling 

of  Monticello  Reservoir  in  South  Carolina.  It  contains  records  of  32 

earthquakes,  recorded  in  May  and  June  of  1979  and  used  in  a  study  of 

stress  drops  by  Fletcher  (1982).  He  reports  they  are  all  magnitude 

M  <  1.7.  Secor  et  aK  (1982)  report  the  earthquakes  as  occurring  in  a 

heterogeneous  quartz  monzonite  pluton  of  Carboniferous  age.  Zoback 

and  Hickman  (1982)  suggest  that  the  earthquakes  are  shallow,  occurring 

along  existing  fractures,  and  the  result  of  an  increase  in  pore 

pressure  caused  by  increasing  the  head  at  the  reservoir.  Figure  1 

adapted  from  Talwani  and  Hutchenson  (1982)  shows  the  location  of  the 

seismic  stations  and  the  reservoirs.  Epicenters  for  the  earthquakes 

considered  have  been  indicated  on  this  figure.  They  are  all  within 

8  km  of  the  recording  stations.  .  Secor  et  al_.  (1982)  report  on  the  two 

reservoirs.  Parr  Shoals  reservoir  was  Impounded  in  1914  and  increased 

in  1976  to  a  volume  of  0.039  km  with  a  surface  area  of  18  knr  and  a 

maximum  depth  of  11  meters.  Monticello  reservoir  has  a  volume  of 
3  2 

0.5  km  ,  a  surface  area  of  27  km  ,  and  a  maximum  depth  of  48  meters. 


Fletcher  (1982)  describes  the  Instrumentation  used  to  record  the 
Monticello  earthquakes.  A  summary  is  given  here.  Sprengnether 
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Table  1.  Locations  of  earthquakes  and  recording  stations  at 
Monticello,  South  Carolina. 


EVENT/STATION  "S"  TIME 


LOCATION 


1271000 

34  20.03  N 

081  19.62  W 

OUC 

0.68 

2.2 

34  20.07  N 

081  21.06  W 

DON 

1.05 

3.5 

34  21.42  N 

081  21.20  W 

LKS 

0.89 

3.0 

34  19.95  N 

081  17.69  W 

1320218 

NOT  POSSIBLE 

JAB 

1.66 

5.8 

34  22.28  N 

081  19.47  W 

1281119 

34  20.70  N 

081  20.81  W 

JAB 

1.31 

4.6 

1302328 

34  20.49  N 

081  20.65  U 

SNK 

0.56 

1.8 

34  20.29  N 

081  19.54  W 

JAB 

1.19 

3.8 

LKS 

1.38 

4.7 

1310603 

34  18.50  N 

081  20.50  W 

SNK 

1.05 

3.7 

LKS 

1.48 

5.2 

JAB 

2.00 

7.0 

1281831 

NOT  POSSIBLE 

JAB 

1.02 

3.57 

- -«  V  >  -•  v  --  -•  V  V  V  V  V  -* 


Table  2.  Locations  of  Earthquakes  and  Recording  Stations  at  Mammoth 
Lakes,  California. 


EVENT/STATION  "S"  TIME 


LOCATION 


1571941 

37  32.91  N 

118  52.53  W 

MGE 

3.17 

8.0 

37  33.67  N 

118  47.22  W 

HCF 

4.37 

10.6 

37  38.51  N 

118  50.98  W 

CBR 

5.50 

15.2 

37  40.75  N 

118  49.51  W 

TWL 

4.84 

13.7 

37  36.93  N 

119  00.37  W 

ROC 

5.67 

14.9 

37  29.78  N 

118  43.16  W 

TOM 

6.11 

18.0 

37  33.05  N 

118  40.32  W 

LKM 

5.97 

17.3 

37  41.80  N 

118  56.13  W 

1592317 

37  37.50  N 

118  52.52  W 

HCF 

2.38 

3.0 

FIS 

2.58 

4.2 

37  36.84  N 

118  49.82  W 

CBR 

3.60 

7.5 

LKM 

3.73 

9.6 

MGE 

4.03 

10.6 

TWL 

4.52 

11.6 

LAK 

5.30 

13.1 

37  38.49  N 

118  43.70  W 

DR-100,  12-bit,  digital  recorders  were  used  with  an  unspecified  make 
of  geophone.  The  sample  rate  was  200  Hz  with  an  antialiasing  filter 
of  either  50  Hz  or  70  Hz  having  a  30  db  per  octave  rolloff.  The 
geophones  were  velocity  sensors  with  damping  specified  as  0.7  of 
critical.  The  possibility  of  aliasing  on  Monticello  records  for  which 
the  filter  comer  frequency  was  set  at  70  Hz  was  considered.  Inspec¬ 
tion  of  bandpass  filtered  records  at  the  higher  frequencies  revealed 
no  difference  in  the  rate  of  amplitude  decay  between  the  recordings 
made  with  a  50  Hz  corner  filter  and  a  70  Hz  corner  filter. 

The  second  set  of  data  is  from  Mammoth  Lakes,  California.  It 
consists  of  the  150  earthquakes  used  by  Archuleta  et  aK  (1982)  in  a 
study  of  source  parameters.  The  earthquakes  followed  a  series  of 
magnitude  M  >  6.0  earthquakes  which  occurred  May  25,  1980  (Spudich  et 
al . ,  1981).  Wallace  et^.  (1982)  Investigate  a  discrepancy  in  fault 
plane  projections  as  reported  in  the  literature  for  the  Mammoth  Lakes 
earthquakes.  Citing  evidence  of  an  active  magma  body,  they  suggest  a 
low  velocity  zone  related  to  recent  magmatic  activity  as  a  possible 
cause.  Figure  2  adapted  from  Spudich  et  al_.  (1981)  shows  the  station 
locations.  Epicenters  for  the  earthquakes  considered  have  been  added 
to  this  figure.  The  source  receiver  distances  range  from  5  to  25  km. 
The  Long  Valley  Caldera,  a  volcanic  crater  0.7  million  years  old,  is 
outlined  along  with  major  faults.  The  Instrumentation  at  Mammoth 
Lakes  Is  essentially  the  same  as  at  Monticello,  except  that  accelero¬ 
meters  are  used  at  some  of  the  Mammoth  stations.  Details  of  the 
instrumental  operations  and  clock  corrections  can  be  found  in  Spudich 
et  a}..  (1981). 


Figure  2*  Recording  Stations  and  Epicenters  at  Mammoth 
Lakes*  California*  stations(A)*  epicenters <*) 


Analysis 


The  rms  coda  amplitude  is  related  to  the  power  spectrum  by 
A(f  |t)  =  [P(f  |t)Af]^.  Substituting  the  equation  for  the  power  spec¬ 
trum  in  the  single  scatter  model  gives 

A(f|t)  =  ct"m  exp(-irft/QT)  .  (9) 

Over  the  range  of  frequencies  for  which  geometric  scatter  is  the  case, 
is  constant.  Therefore,  c  =  (sAf)*^  has  the  same  frequency 
dependence  as  the  source  for  octave  Af.  For  cylindrical  waves 
m  =  0.5,  while  for  spherical  waves  m  =  1.0.  Taking  the  ordinary 
logarithm  of  this  expression 

log[A(fjt)]  =  C  -  m  .  log(t)  -  bt  (10) 

is  obtained  where  b  *  log(e)irf/QT  and  C  *  log(c) .  To  estimate  QT  this 
equation  is  fit  to  filtered  records  by  the  method  of  least  squares. 
This  is  the  approach  taken  by  Aki  and  Chouet  (1975) . 

To  study  the  temporal  decay  of  coda  waves  the  longest  records  are 
chosen.  Of  these  the  duration  of  recording  ranges  from  8.0  to  23.0 
seconds  at  Monticello  and  from  10.0  to  35.0  seconds  at  Mammoth  Lakes. 
Only  a  few  records  from  Monticello  are  considered  to  be  of  sufficient 
duration  for  coda  analysis  (Table  1).  The  Mammoth  Lakes  records 
chosen  for  analysis  are  listed  in  Table  2. 

Travel  times  are  taken  from  the  location  data  provided  with  the 
U.S.G.S.  tapes.  Some  of  the  Monticello  earthquakes  were  not  located. 
For  these  the  S-wave  travel  time  is  calculated  from  the  difference  in 
S  and  P  wave  arrival  times  and  assumed  S  and  P  wave  velocities  of 


3.5  km/s  and  6.0  km/s  respectively.  Ground  displacement  amplitude  as 
a  function  of  time  in  octave  wide  bands  is  determined  for  the  selected 
records  by  the  following  procedure. 

The  mean  is  removed  from  a  record  in  a  256  sample  (1.28  second) 

window  centered  about  the  origin  time.  This  portion  of  the  trace  is 

multiplied  by  normalized  Hamming  coefficients.  The  windowed  trace  is 

corrected  for  instrument  gain,  Fourier  transformed,  and  multiplied  by 
9 

l/f  or  1/f  .  Assuming  a  flat  instrument  response  in  the  frequency 
band  of  interest,  this  converts  the  trace  to  ground  displacement 
amplitude  from  velocity  and  acceleration  records  respectively.  It  has 
been  assumed,  as  in  Archul eta  et  al_.  (1982)  and  Fletcher  (1982),  that 
the  velocity  response  is  flat  from  2  Hz  to  50  Hz  and  the  acceleration 
response  is  flat  above  0.1  Hz.  The  average  amplitude  is  calculated 
within  octaves  (equal  power  bands)  about  selected  center  frequencies. 
The  log  of  this  amplitude  is  plotted  at  time  zero  on  separate  plots 
for  each  frequency.  The  window  is  then  moved  100  samples  (0.5 
seconds)  and  the  process  is  repeated  with  the  log  average  amplitude 
plotted  at  time  0.5.  This  is  done  until  the  end  of  the  record  is 
reached.  The  result  is  a  series  of  plots  of  log  average  amplitude 
versus  time,  measured  from  origin  time,  at  selected  center  frequencies 
from  3  Hz  to  72  Hz.  Figure  3  contains  a  sample  record  and  its  corres¬ 
ponding  filtered  amplitude  plots.  Frequencies  above  the  corner  of  the 
antialiasing  filter  are  considered  as  the  constant  reduction  in  ampli¬ 
tude  will  only  show  up  in  the  source  term.  Therefore  Qy  can  still  be 
estimated  provided  there  is  sufficient  energy  present  at  these 
frequencies. 
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Figure  3  ,  b.(top).  c. (bottom) 
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Figure  3  .  d.(cop)  e. (bottom) 


00  1  :  3  4  !  4  7  «  9  O  i  2  1  n  4  ?  M  0  1  2  n  5  4  7  I  »  0  1  !  1  «  J  t  7  I 
0  ♦ 

0 


0 

0 

0 

0 

0 

0 

-3. SO  0 


0 

0 

0 

0 

0 

0 

0 

0 

-4.30  0 


♦♦  ♦  *♦♦♦ 


00  123434789012343478901 
LAPS*  TIHC  ( SECONDS  > 

34347840  1  2345478901  23 


34347840123434784 


00  1 
0 


3  4  7  8  4  0  1 
♦ 


3  4  3  4  7  8  4 


-3.80  0 
0 
0 
0 


0 

0 

0 

0 

0 

-4.80  O 
0 
0 
0 
0 


►  ♦  ♦♦♦♦  f 


♦  ♦♦♦♦  ♦♦♦ 


00123494784012343478401234347890 
LAPS*  TIME  (SECONDS) 


123434789 


Figure  3-  h.(top)  i. (bottom) 


The  next  step  is  to  fit  the  model  equation  for  log[A(f|t)]  to  the 
backscatter  portion  of  the  logCA(fjt)]  plots  for  an  estimate  of  QT. 
The  method  of  least  squares  is  used  and  two  fits  are  made  corres¬ 
ponding  to  the  two  choices  of  m  =  1.0  and  m  =  0.5.  Standard  devia¬ 
tions  are  also  calculated  for  each  fit.  Attempts  to  fit  the  model 
taking  m  as  an  independent  variable  produced  erratic  results  as  has 
been  reported  in  the  literature  (Aki  and  Chouet,  1975). 

The  backscatter  portion  of  the  filtered  data  is  determined  as 
follows.  Energy  arriving  after  times  the  primary  S-wave  travel 
time  is  backscatter  in  the  sense  that  the  scattering  angle  is  greater 
than  90  degrees.  However,  because  of  the  1.28  second  window  and  0.5 
second  steps  used  for  averaging  only  arrivals  1.0  second  later  than 
this  are  entirely  backscatter. 

On  some  records  the  dynamic  range  of  the  recording  instrument  has 
been  exceeded  for  the  early  coda.  The  saturated  portion  of  the  coda 
is  not  considered  in  the  analysis.  Also  excluded  are  the  last  three 
points  of  the  filter  output  as  they  are  distorted  by  artificial  exten¬ 
sion  of  the  end  of  the  trace  necessary  In  filtering. 

Another  consideration  Is  noise.  The  background  noise  level  is 
assumed  constant  and  taken  equal  to  the  peak  noise  level  recorded 
prior  to  the  P-wave  arrival.  If  no  noise  is  recorded  here  it  is 
assumed  to  be  +  or  -  0.5  times  one  unit  of  digitization.  The  noise 
level  Is  converted  to  ground  displacement  amplitude.  A  reference,  five 
times  (14db  above)  this  noise  level.  Is  taken  as  the  minimum  accept¬ 
able  signal  level.  Portions  of  the  coda  for  which  the  average  signal 
amplitude  falls  below  this  reference  are  not  considered  in  the  fits. 
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Finally,  from  the  longest  records  it  Is  evident  that  Q-j.  changes 
with  time  causing  the  Q^.  value  obtained  to  depend  on  the  portion  of 
coda  fit.  This  can  be  seen  on  the  sample  record  in  Figure  3.  This 
time  dependence  of  Qj  was  also  observed  by  Raytlan  and  Khalturln 
(1978)  who  suggest  variations  of  Q  with  depth  and  temporal  variations 
of  the  geometric  spreading  factor  as  possible  causes.  The  variation 
in  spreading  factor  could  result  from  a  change  in  the  type  of  waves 
being  recorded  as  a  function  of  time  (e.g.,  surface  to  body  waves  as 
time  increases).  Another  possibility  is  that  of  a  significant  deter¬ 
ministic  influence  on  the  energy  arriving  at  the  receiver.  Long  and 
Wilson  (1982)  have  noticed  resonance  peaks  in  the  power  spectra  of  the 
Monticello  records.  The  frequencies  at  which  these  occur  are  in  con¬ 
gruence  with  conceivably  supportable  resonances  in  the  Monticello 
reservoir.  This  can  be  regarded  as  an  example  of  multiple  scattering, 
which  tends  to  raise  the  value  of  Q  because  energy  previously  scat¬ 
tered  out  of  the  coda  is  scattered  back  in,  reducing  the  rate  of 
decay.  Indeed,  in  the  limit  of  strong  scattering,  there  is  no  contri¬ 
bution  to  Q  from  scattering  (Dainty  and  Toksoz,  1977).  Because  of  the 
observed  temporal  dependence  of  Qj,  fits  are  made  to  different 
portions  of  the  coda  where  possible. 

Results 

Plots  of  the  average  over  all  stations  and  events  of  1  ogCQjjl 
versus  log[frequency]  are  constructed  and  appear  as  Figure  4.  Because 
a  temporal  dependence  of  Qj  was  observed,  three  plots  are  made  for 
each  area  and  m  value.  They  are  from  the  beginning  of  backscatter  to: 
I,  five  seconds  beyond:  II  between  five  and  nine  seconds  beyond;  III, 
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more  than  nine  seconds  beyond.  These  times  were  chosen  arbitrarily 
within  the  limits  of  the  data.  As  the  plots  are  roughly  linear,  best 
fits  of  the  equation  Qj  *  afb  are  made  to  each  of  these  to  determine 
the  values  of  a  and  b.  This  approach  Is  common  In  the  literature. 
The  a  and  b  values  obtained  are  presented  In  Table  3. 

The  frequency  dependence  of  Qy  may  be  used  to  estimate  Q^  and  g 
If  some  assumptions  are  made  about  their  frequency  dependence.  For 
the  single  backscatter  model  applied  here  the  density  of  scatterers  p 
Is  assumed  constant.  The  backscatter  cross  section  ob  will  be  a  con¬ 
stant  for  geometric  (high  frequency)  scatter.  The  result  is  a  con¬ 
stant  g  for  geometric  scatter  with  the  single  backscatter  model.  It 
Is  usually  assumed  that  Q^  Is  Independent  of  frequency  in  agreement 
with  most  laboratory  measurements  of  (Knopoff,  1964).  With  these 
assumptions  a  least  squares  fit  of  1/Qy  *  1/Q^  +  gB/Zwf  Is  made  to  the 
Qy  versus  frequency  data  to  estimate  g  and  Q1  (Dainty,  1981).  Nega¬ 
tive  values  of  occur  when  the  least  square  fit  results  in  a  gB/2wf 
which  is  greater  than  1/Qy.  As  Increases  with  Increasing  gB/2irf  a 
negative  result  Is  Indicative  of  a  very  high  value.  These 
results  are  shown  In  Table  4. 

A  second  estimate  of  g  is  obtained  as  follows.  From  the  least 
squares  fit  of  log[A(f  t)]  estimates  of  Qy  and  the  source  factor  have 
been  obtained.  The  filter  output  at  the  primary  S-wave  travel  time  at 
a  given  frequency  Is  a  component  of  the  source  factor  C.  This  allows 
a  second  estimate  of  g  from  the  primary  S-wave  spectra  and  the  source 
factor  C  ■  log  (saf). 

Taking  a  result  from  the  theory  section 


Table  4 

I  _ 


Log(g)  Values  From  Fitting  Equation  (5). 


Monticello 


Mammoth 


51 


s  =  |*(f|rQ)|2  4mirr02m  (B/2)1’2m  g  .  (11) 

Also,  the  average  amplitude,  |*(f|r)|,  of  waves  incident  at  a  scat- 


terer  at  distance  r  is 


| ♦  ( f  |r)  |  =  [♦(f  |rQ)|  (r/rQ)"m  exp(-Trfr/BQT)  .  (12) 

It  is  important  to  note  that  |$  (f|r)|  is  also  the  average  spectral 
amplitude  of  the  direct  S-wave  where  r  is  the  source  receiver  dis¬ 
tance.  This  is  estimated  from  a  linear  interpolation  of  the  spectral 
amplitudes  measured  in  the  two  windows  closest  to  the  S  arrival  time. 
The  interpolation  is  necessary  as  the  window  is  generally  not  centered 
on  the  S  arrival.  Solving  equation  (12)  for  the  source  spectrum  at 
r0,  j *  (f|r0)|,  and  substituting  in  (11)  a  solution  for  log(g)  is 
obtained. 

log(g)  *  2C  -  log[4rnir(B/2)1"2m  j  ♦  (f j r) j 2  r2m  exp(2irft/Qy)Af  ]  . 


Estimates  of  average  log(g)  are  obtained  using  this  equation  for 
m  =  0.5  and  or  m  *  1  and  are  tabulated  in  Table  5. 


Discussion 

The  log(Qy)  versus  log (frequency)  curves  of  Figure  4  show  Qy  to 
be  roughly  the  same  for  similar  fits  at  Mammoth  and  Monticello. 
Though  it  is  generally  agreed  that  Qy  Is  higher  in  the  eastern  U.S. 
the  results  here  may  be  explained  by  the  short  lapse  times  of  the 
codas  used  to  estimate  Qy,  These  lapse  times  are  much  shorter  than 
are  usually  reported  in  the  literature  for  Qy  from  coda  estimates. 


Table  5.  Log(g)  as  Estimated  from  the  S-wave  Spectrum  and  Coda  Source 
Function. 


This  Implies  that  the  results  here  represent  the  near  surface  Qy  and 
that  It  Is  approximately  the  same  for  both  regions. 

The  results  presented  In  Table  3  are  the  a  and  b  values  obtained 
by  best  fit  of  Qy  »  af*5  to  the  log[Qy]  versus  log[ frequency]  curves  of 
Figure  5.  The  values  of  b  In  Table  3  range  from  0.6  to  1.8,  with  an 
average  of  1.3.  The  usual  value  of  b  reported  In  the  literature  Is 
around  0.5  (0.6,  Fedotov  and  Boldyrev,  1969;  0.5,  Rautlan  and 
Khalturln,  1978;  0.5,  Tsujiura,  1978;  0.6,  Akl,  1980).  However,  the 
higher  values  obtained  here  are  not  unique.  Rovelll  (1982)  has 
reported  a  b  value  of  1.1  for  Friuli,  Italy.  This  value  was  determined 
by  coda  analysis.  Akl  suggests  that  the  b  value,  0.8,  which  he  deter¬ 
mined  for  northeast  Kanto,  Japan,  Is  the  result  of  Increased  scat¬ 
tering  by  a  higher  density  of  scatterers  associated  with  the  higher 
level  of  tectonic  activity  In  this  part  of  the  region.  Aki  and  Chouet 
(1975)  report  on  Q  versus  frequency  at  Stone  Canyon,  California.  A  b 
value  of  0.9  is  obtained  by  best  fit  of  Qy  =  afb  to  the  reported 
results.  Stone  Canyon,  Mammoth  Lakes,  and  Friuli  are  all  highly 
active  regions. 

The  b  value  results  at  Montlcello  may  be  considered  anomalous  for 
this  region  of  relative  tectonic  stability.  However,  in  addition  to 
Including  the  reservoirs  as  a  substantial  part  of  the  region  sampled 
the  codas  and  hypocentral  distances  at  Montlcello  are  much  shorter 
than  are  usually  reported  In  the  literature  for  other  studies.  Conse¬ 
quently  the  near  surface  scatterers  are  favored  in  sampling.  Andrews 
(1982)  has  suggested  that  near  surface  scattering  is  the  dominant 
scattering  effect  for  the  similar  situation  at  Mammoth  Lakes. 


Table  4  gives  the  results  of  fitting  equation  (5)  to  the  values 
of  Q-j.  as  a  function  of  frequency^  Qi  and  g  are  considered  constants 
with  frequency  In  making  this  fit.  In  (5),  the  two  terms  on  the  right 
hand  side  will  be  equal  for  a  frequency  F  given  by 

F  -  gBQ1/(2w)  (14) 

Table  6  gives  the  values  of  F  found  using  the  results  presented  In 
Table  4.  For  frequencies  less  than  F,  scattering  dominates  the 
attenuation,  while  for  frequencies  greater  than  F,  inelastic 
attenuation  will  dominate.  From  Table  6,  scattering  dominates 
attenuation  over  most  of  the  frequency  range  in  most  cases,  especially 
at  Montlcello.  From  (5),  if  scattering  dominates  then  QT  should  be 
proportional  to  frequency,  In  approximate  agreement  with  the  results 
presented  In  Table  3.  The  proportionality  of  QT  with  frequency  Is 
also  seen  In  Figure  4. 

There  Is,  however,  another  issue  raised  by  Table  6.  Dainty  and 
Toksoz  (1981)  propose  that  If  the  frequency  f  <  F,  multiple  scattering 
Is  a  possibility.  Obviously  this  situation  holds  for  most  of  the 
codas  considered  here,  from  Table  6.  Dainty  and  Toksoz  point  out  that 
this  seems  to  be  true  for  many  cases  cited  In  the  literature. 
However,  Dainty  and  Toksoz' s  condition  is  only  a  necessary,  not  a 
sufficient,  condition  for  multiple  scattering.  Since  It  will  take 
some  time  after  the  origin  time  for  multiple  scattering  to  become 
Important,  a  condition  from  Sato  (1977)  that  multiple  scattering  will 
become  Important  for  times  t  >  T,  where  T  *  l/(gB),  Is  appropriate;  If 
f  <  F  and  t  >  T,  strong  scattering  Is  probably  occurring.  T  is  also 
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Table  6.  F  =  960,7(211),  T  =  l/(gB).  *  is  F  tending  to  infinity. 


Monti cello  Mammoth 


F 

Hz 


30 

15 

20 
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given  In  Table  6,  and  In  general  Indicates  that  multiple  scattering 


might  be  occurring  in  some  of  the  codas,  remembering  that  I  Indicates 
t  _<  5  sec;  II,  t  9  sec;  III,  t  >  9  sec;  but  In  practice  t  Is  less  30 
sec  for  all  codas  In  III.  More  generally,  both  F  and  T  In  Table  6  are 
of  the  same  order  of  magnitude  as  the  frequencies  and  times  considered 
In  this  study,  suggesting  that  the  situation  may  be  transitional 
between  weak  and  strong  scattering.  A  similar  conclusion  was  reached 
by  Dainty  and  Toksbz  (1981).  Theoretical,  formulations  to  handle  this 
problem  (e.g.,  Gao  et  al . ,  1983)  do  not  appear  to  conserve  energy  and 
do  not  appear  to  agree  with  the  strong  scattering  diffusion  formalism 
of  Dainty  and  Toksbz  (1977). 

The  log[g]  values  obtained  from  the  coda  source  spectrum  are 
presented  In  Table  5.  In  general  they  are  about  an  order  of  magnitude 
to  two  orders  of  magnitude  smaller  than  the  log[g]  values  obtained 
from  Dainty's  equation.  Andrews  (1982)  has  found  the  source  spectral 
amplitude  from  the  coda  to  be  ten  times  larger  than  the  spectral 
amplitude  measured  from  the  S-wave.  As  the  value  of  g  here  is 
proportional  to  the  quotient  (coda  source  spectrum/S  wave  source 
spectrum),  a  small  g  value  means  our  results  do  not  agree  with  those 
of  Andrews'.  However,  they  are  In  agreement  with  the  theoretical 
results  reported  In  the  second  part  of  this  report. 

Throughout  this  report,  the  question  of  whether  the  scattering  is 
two-dimensional  (m  *  0.5,  cylindrical  spreading,  surface  waves)  or 
three-dimensional  (m  *  1.0,  spherical  spreading,  body  waves)  has  been 
left  open.  As  discussed  in  the  section  on  Analysis,  It  was  not 
possible  to  determine  the  correct  value  of  m  from  the  data.  Akl 
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(1980)  has  suggested  on  various  grounds  that  the  three-dimensional 
case  Is  probably  appropriate,  but  he  dealt  with  much  longer  codas. 
Most  of  the  more  Important  results  presented  here  seem  to  hold  for 
both  cases.  It  must  be  conceded,  however,  that  the  change  in  coda 
decay  as  time  Increases  could  be  due  to  a  change  in  QT  reflecting  the 
change  from  surface  waves  to  body  waves. 
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High-Frequency  Acoustic  Backscattering  and  Seismic  Attenuation 

Anton  M.  Dainty 

School  of  Geophysical  Sciences,  Georgia  Institute  of  Technology 

An  expression  lor  the  backscattered  intensity  of  acoustic  waves  singly  scattered  from  a  region  contain¬ 
ing  fluctuations  of  the  acoustic  velocity  has  been  derived  for  high  frequencies  by  expanding  the  auto¬ 
correlation  of  the  slowness  fluctuations  in  a  Taylor  series  about  aero  lag.  The  resulting  expression 
indicates  that  the  backseat  tend  intensity  is  independent  of  frequency  and  directly  proportional  to  the 
first  derivative  of  the  autocorrelation  at  zero  lag;  the  next  higher  term  is  proportional  to  the  reciprocal  of 
the  square  of  the  frequency  and  directly  proportional  to  the  third  derivative  of  the  autocorrelation  at 
aero  lag.  Contributions  from  terms  of  the  Taylor  series  involving  even  numbered  derivatives  of  the 
autocorrelation  sre  zero.  Since  for  the  autocorrelation  of  s  smooth  function  only  the  even  numbered 
derivatives  sre  nonzero  at  zero  lag,  this  result  demonstrates  that  backscattering  at  high  frequencies  can 
only  oocur  from  discontinuities  of  velocity  or  its  derivatives  u  opposed  to  fluctuations  in  which  the 
velocities  are  smooth.  If  the  backscattered  intensity  is  independent  of  frequency,  the  contribution  of 
backscattering  to  the  attenuation  parameter  Q  is  proportional  to  frequency.  Such  behavior  may  have 
been  observed  for  seismic  waves  of  frequencies  greater  than  1  Hz,  suggesting  that  scattering  from 
discontinuities  is  an  important  part  of  the  attenuation  of  such  waves. 


Introduction 


Random  Medium  Models 


The  phenomenon  of  scattering  of  seismic  waves  has  attrac¬ 
ted  increasing  attention  in  recent  years.  Backscattering  of  seis¬ 
mic  waves  has  been  used  to  explain  terrestrial  codas  [Aid  end 
Ckouet,  1973;  Dainty  and  Toksot,  1981]  and  part  or  all  of  the 
attenuation  suffered  by  seismic  waves  at  frequencies  of  1  Hz 
and  greater  [Aki,  1980;  Dainty,  1981;  Kikuchi,  1981;  Wu, 
1982].  To  analyze  observed  data  on  attenuation,  coda  ampli¬ 
tudes  and  other  manifestations  of  scattering,  a  model  of  the 
earth  as  a  scattering  medium  is  required  together  with  a 
theory  to  derive  quantitative  results  from  this  model.  Two 
general  types  of  model  have  been  used.  One  is  the  “randomly 
distributed  specific  scatterers”  model  This  model  assumes 
scatterers  of  a  specific  type,  for  example,  spheres,  randomly 
distributed  throughout  an  otherwise  homogeneous  region.  Ex¬ 
amples  of  this  type  at  model  are  the  sphere  model  of  Dainty 
[1981],  used  to  explain  the  frequency  dependence  of  the  at¬ 
tenuation  parameter  Q  at  frequencies  above  1  Hz,  and  the 
crack  model  of  Kikuchi  [1981],  also  used  to  analyze  Q.  A 
second  type  of  model  is  the  “random  medium"  model  in 
which  the  earth  is  considered  to  have  average  properties  that 
do  not  vary  from  place  to  place  but  with  random  fluctuations 
of  the  material  properties.  Examples  are  the  use  of  an  acoustic 
random  medium  by  Aid  [1980]  and  Wu.  [1982]  to  analyze  Q 
above  1  Hz. 

In  this  paper  I  shall  examine  the  acoustic  random  medium 
model  using  the  results  of  Chernov  [I960]  in  the  limit  of  high- 
frequency  backscattering.  This  problem  was  chosen  because 
investigators  using  this  model  have  generally  concluded  that 
this  is  the  practical  case  in  the  earth  for  frequencies  of  1  Hz 
and  higher  [Aid.  1980;  Dainty,  1981;  Sato,  1982 b\  Wu,  1982]. 
By  analyzing  the  high-frequency  limit,  insight  may  be  gained 
as  to  what  type  of  information  about  the  medium  may  be 
deduced  from  high-frequency  data.  An  additional  benefit  is  an 
understanding  of  the  relationship  between  the  two  types  of 
models  that  have  been  used  to  examine  this  problem. 
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The  simplest  type  of  randQja-Tfiedium  used  in  seismology  is 
the  acoustic  random  medium  discussed  by  Chernov  [I960]  (a 
summary  of  some  useful  results  is  given  by  Aid  and  Richards 
[1980]).  The  medium  is  considered  to  have  fluctuations  of 
velocity  and  density  about  mean  values;  for  simplicity  I  will 
consider  only  the  velocity  fluctuations.  Let  the  acoustic  veloci¬ 
ty  C  be  given  by 

C-C„[l-M*)]  (1) 

where  C0  is  a  constant  and  the  slowness  perturbation,  is 
a  random  function  of  position  x  and  has  an  average  value  of 
zero.  The  autocorrelation  N(t)  of  the  slowness  perturbation  it 
is  then 

N{r)  «  <Mx)/4*  +  r»/<y>  (2) 

The  vector  r  is  the  lag,  and  the  angle  bracken  indicate  averag¬ 
ing  over  all  x.  If  the  fluctuation  is  isotropic,  then 

Nit)  -  jV(r)  (3) 

that  is,  the  autocorrelation  depends  only  on  the  magnitude  of 
the  lag  and  not  its  direction. 

Chernov  [1960]  finds  the  amplitude  of  singly  (weakly)  scat¬ 
tered  waves  for  the  case  of  a  plane  wave  traveling  through  an 
acoustic  medium  of  velocity  C0  without  velocity  fluctuations, 
except  for  a  volume  V  containing  velocity  fluctuations  as  de¬ 
scribed  above.  If  the  incident  (pressure)  wave  is  given  by 


Pi  -  A  •  exp  [2x/?(f  -  x/C„)]  (4) 


where  /  is  the  frequency,  then  the  amplitude  of  the  scattered 
wave  at  position  R  from  the  center  of  V  in  the  limit  of  R  large 
in  comparison  with  the  linear  dimension  of  V  is  given  by 


4 

Co5*1  sin  (9/2) 


■  sin  {[4x/r  ■  sin  (0/2)]/Co}  rdr  (5) 


{Chernov,  1960,  equation  43,  p.  31].  The  scattering  angle  9  is 
the  angle  between  R  and  the  direction  of  propagation;  I  shall 
define  “backscattering”  as  the  case  9  >  x/2.  Equation  (5)  was 
first  derived,  in  a  slightly  different  form,  by  Pekeris  [1947]. 

3172 


•V.  ‘  1-  A  - 


Dainty:  Hjoh-Fuqusncy  Acoustic  Backscattbuno 


This  paper  is  concerned  with  the  evaluation  of  (5)  in  the  limit 
of  large/and  0  >  a/2. 

1  Before  proceeding,  it  should  be  noted  that  the  acoustic 
medium  described  above  is  obviously  too  simple  to  describe 
the  earth,  which  is  an  elastic  medium  presumably  containing 
fluctuations  of  both  elastic  constants  and  density  and  through 
which  two  types  of  waves,  compressions!  and  shear,  can  prop¬ 
agate.  Knopoff  and  Hudson  [1967]  have  demonstrated  that 
converted  waves  (scattered  compressions!  waves  from  incident 
shear  waves,  and  vice  versa)  are  not  important  at  high  fre¬ 
quencies.  Haddon  and  Cleary  [1974]  examined  the  problem  of 
scattering  from  an  inhomogeneous  elastic  medium;  some  of 
their  (unpublished)  results  are  discussed  in  the  work  of  Akt  and 
Richards  [1980].  These  results  indicate  that  at  least  some  of 
the  formulas  equivalent  to  equation  (5)  are  of  a  similar  form. 
Accordingly.  (5)  may  be  considered  typical  of  integrals  that 
appear  in  this  type  of  problem. 

Evaluation  of  the  Backscattered  Intensity 
at  High  Frequency 

To  evaluate  (5),  we  need  to  find  the  integral 

f" N{r)  ■  sin  {[Anfr-an  <9/2)]/C0}  ■  rdr  -  F(ao)  -  F( 0)  (6) 

where  F  is  the  indefinite  integral  However,  if  there  is  a  scale 
length  a  associated  with  the  fluctuations,  N(r)  will  generally 
decline  rapidly  for  r  larger  than  a.  Accordingly,  I  will  assume 
that  F(oo)— »0  in  (6);  this  will  certainly  be  true  of  Af(r)-» 
exp  [— r/a]  for  large  r.  Then 

N(r)  -  sin  {[4 njr  -  sin  (0/2)]/C*}  rdr  a  -  F(0)  (7) 

To  approximate  F(0)  in  the  cass  of/large  and  0  >  x/2,  note 
that  for  this  case  the  sine  term  in  (7)  varies  much  more  rapidly 
than  N{r\  assuming  N  is  a  smooth  function.  Accordingly,  I 
will  expand  N(r)  in  a  Taylor  series  near  r  -  0.  First,  however, 
define 

«-[4*frsin(fl/2)]/C0  (8) 


j* u-  sin  (u)-du  =*  — u"  cos  (u)  +  m  j* u"'1  cos  (u)-du  (14) 

j"  u*  cos  (u)  •  du  m  um  sin  (u)  —  m  j* u"~ 1  sin  (u)  du  (15) 

J* u- cos  (u)du  m  cos  (u)  +  u- sin  (u)  (16) 


•cos  (u)  du  —  cos  (u)  +  u  sin  (u) 


Using  (14)  and  (16), 


sin  (u)  <fu  »  2u  sin  (u)  +  (2  -  uJ)  cos  (u)  (17) 


Using  (13).  (14),  (15),  and  (16), 


r  r  c»-iwi 

I  sin  (u)  du  —  cos  (i»)  |  —  u"  +•  £  (  — If*1 

•  m(m  —  1)  •  ■  •  (m  —  2n  +  l)u*_2*J 

r  lm-U/2 

+  sin  (u)|  «u“* 1  +  £  (-1)" 

m(*-  l)...(m-2n)u— 


(m  odd) 


j*iT  sin  (u)-du  —  cos  (u£— u"  +  £(  —  l)**' 

m(m  —  1) ...  (m  —  2n  + 

r  m-i’ 

+  sin  (us  nuT~ 1  +  £  (-1)* 

L  —i 

l)...(m-2»)ir-2*-,J 


(m  even) 


*C  -  [4*/-  sin  (0/2)]/Co 


Then  (7) 


From  (131(17),  (18),  and  (19).  as  u-»0. 


p  J[" *(£)  ‘  “  »•«*•*  ~f(°>  U<* 

Expand  H{u/K)  as,  remembering  that  u/K  -  r, 

«(i) (f)^  ♦  i  ( j)  (£?),.„  *  •  •  •  (1» 

provided  the  derivatives  at  r  ■  0  exist  and  are  finite.  Then 

w“p[wJ’-*(**+z(SL 

•  Ju*  sin  (u)'du  +  "J  (12) 

To  evaluate  (12),  use  the  tabulated  integrals  [Comdchoef  and 
Smith,  1962]: 

fu-sin(u)'duwsin(H)-u-cos(u)  (13) 


>  0  (m  odd) 


Jusin  («)dn-»0  (20) 

Ju2  sin  (u)-du-* 2  (21) 

J*u"  sin  (u)-du->0  (modd)  (22) 

Ju"  sin  (u)  du-»(-tr'2*',m!  (m  even)  (23) 
Substituting  in  (12)  and  letting  u—  0, 

„0).  y  ‘  ,.r-W  a  (24) 

W  r  (2n  -  1)!  V  dr»- 1  )r.0 

According  to  (24),  under  the  conditions  described  the  back- 
scattered  intensity  depends  only  on  the  odd  derivatives  of  N  at 
r  -  0  and  not  on  the  even  derivatives.  Since  this  is  a  high- 
frequency  approximation,  the  first  term  in  (24)  will  be  the 
most  important,  that  is. 


(dN\ 

'\  +  Lt 


+  0d/K») 
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Substituting  in  (10)  and  (3), 

M2</i2> 


w*  -  - 


_ « 

8it R1  sin*  (9/2)  \  dr  ), 


■0(1  /f1)  (26) 


This  indicates  that  the  backscattered  intensity  is  independent 
of  frequency,  provided  that  the  autocorrelation  has  a  nonzero 
derivative  at  zero  lag.  Since  the  first  derivative  at  zero  lag 
must  be  negative  (N(r)  has  a  maximum  at  r  »  0),  |p,|2  will  be 
positive. 

To  test  (26),  I  will  compare  it  to  two  specified  cases  of  JV(r) 
for  which  explicit  formulas  for  |p,j2  are  available  [ Chernov , 
I960],  If 

N(r)  «  exp  (-r/a)  (27) 

then 


32  nW/VQi2) 

W  ~  R2[C02  +  16ir2/2a2  sin2  (6/2)] 2 


(28) 


\P,Y 


8 xR2a  ■  sin*  (6/2) 


f— *  oo,  6  2  */2 


(29) 


The  approximate  solution  given  by  (26)  reproduces  this  result 
Another  case  for  which  an  exact  solution  is  available  is 


(V(r)-exp(- r2/a2)  (30) 

Then 

2  4xV2VA1f*a1(ft2) 

W  "  R2 

exp  {-[4*2/2a2  sin1  <8/2)]/C02}  (31) 

IpJ  2  —  0  f  ►  oo  6  ^  n/2)  (32) 

Since  N(r)  has  no  nonzero  derivatives  of  odd  order  for  r  =  0, 
this  is  also  in  accord  with  (26). 


Discussion 

The  limit  of  validity  of  (26)  can  be  obtained  from  (7),  since 
we  require  that  the  sine  terms  in  (7)  vary  much  more  rapidly 
than  N(r).  If  a  is  a  scale  length  describing  the  range  of  r  over 
which  N(r)  varies  appreciably,  we  must  have 

[4*/a  •  sin  (6/2 )]/C„  »  2a  (33) 

or 

[2/a  •  sin  (8/2)]/C0  »  1  (34) 

Note  that  (34)  requires  that  the  scattering  angle  6  cannot  be 
small;  in  fact, 

sin  (8/2) »  C0/(2fa)  9  a  30°  (35) 

is  necessary.  Thus  the  approximation  used  here  is  not  valid  for 
forward  scattering.  It  may,  however,  be  valid  for  8  <  a/2  pro¬ 
vided  the  restriction  (35)  is  observed,  together  with  a  re¬ 
striction  discussed  below  requiring  6  s  30°.  Within  these 
limits,  the  scattered  intensity  varies  as  l/(sin*  6/2),  indicating 
that  “side  scattering"  (8  ~  a/2)  is  about  4  times  as  strong  as 
“backscattering”  (8  -  a)  in  the  strict  sense. 

There  is,  however,  another  consideration  that  must  be  ap¬ 
plied  to  the  forward  scattering  question.  In  forward  scattering 
(8-*0)  the  scattered  wave  merges  with  the  original  incident 
wave  [Chemoe,  1960;  Knopoff  and  Hudson,  1967],  The  theory 
usfd  to  derive  (5)  does  not  satisfactorily  account  for  this 
IChtmot,  1960]  and  accordingly  (26)  should  never  be  used  for 


small  6,  (33)  notwithstanding.  Sato  [1982 a],  in  fact,  indicates 
that  the  merging  of  the  scattered  and  incident  wave  is  a  signifi¬ 
cant  effect  for  6  <,  30°  (see  (44)  and  (45),  below),  and  thus  I 
have  added  the  restriction  8  2  30°  to  (35).  It  should  also  be 
noted,  however,  that  if  ( dN/dr)„0  *  0,  (26)  indicates  that  there 
is  always  backscattered  energy,  no  matter  how  high  the  fre¬ 
quency  of  the  incident  wave.  I  shall  examine  the  type  of  veloci¬ 
ty  fluctuation  that  gives  rise  to  such  a  situation  below. 

An  important  application  of  the  theory  of  backscattering  is 
the  interpretation  of  the  frequency  variation  of  Q  for  shear 
waves  in  the  range  1-30  Hz  [<4fci,  1980;  Wu,  1982].  Q,  the 
quality  factor,  describes  the  attenuation  of  the  amplitude  A  in 
(4)  by 

A(x)  -  A(0)  exp  [— ir /x/(gC0)]  (36) 

The  total  Q  represents  the  effects  of  energy  loss  from  the  pri¬ 
mary  wave  by  the  combined  mechanisms  of  loss  of  seismic 
energy  to  heat  (&)  and  loss  due  to  backscattering  (Qs)  accord¬ 
ing  to  the  relation  [Warren,  1972] 

i/e  -  i/e.  + 1 iQ.  07) 


Wu  [1982]  derives  an  expression  for  Q„  with  the  proper 
choice  of  V,  as 


1/0,  -  Jje  ipj2  sin  (6)  (38) 

where  d>  is  an  azimuthal  angle  about  the  direction  of  propaga¬ 
tion  of  the  primary  wave  and  the  integral  over  scattering  angle 
6  is  taken  over  all  backscattering  angles. 

Substituting  for  IpJ2  from  (26)  and  performing  the  integral 
over  </>,  in  the  limit  of  high  frequency. 


I/O  C0<p2>  fdN\  p  sin  (6)  d9 
8 nf  \dr),m0  J.,2  sin*  (8/2) 

This  may  be  evaluated  using  the  trigonometric  identity 

2  sin2  (8/2)  —  1  —  cos  (8) 


(39) 

(40) 


Comparing  this  with  an  expression  from  Dainty  [1981], 


t IQ,  »  gC0/(2*f)  (42) 

then 


This  indicates  that  g,  the  turbidity,  is  independent  of  frequency 
and  depends  on  the  product  of  the  mean  square  slowness 
fluctuation  and  the  derivative  of  the  autocorrelation  of  the 
slowness  fluctuation.  The  lack  of  dependence  of  g  on  fre¬ 
quency  is  the  same  result  obtained  by  Dainty  [1981]  for  a  set 
of  randomly  distributed  spheres  in  geometric  scatter;  thus  we 
may  say  that  the  approximation  explored  in  this  paper  is  the 
analog  of  geometric  scattering  from  obstacles.  Indeed,  obser- 
vationally  we  cannot  distinguish  between  these  cases,  nor  can 
we  distinguish  between  different  autocorrelation  functions  N(r) 
except  insofar  as  their  derivatives  at  r  -  0  differ.  The  lack  of 
dependence  of  g  on  frequency  may  be  experimentally  observed 
for  seismic  waves  above  1  Hz  [Dainty,  1981]. 

In  (38)  the  integral  over  the  scattering  angle  was  taken  over 
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all  backscattered  angles.  Sato  [1982a]  has  used  a  different 
approach  (the  mean  wave  method)  from  that  used  by  Chernov 
[1960]  and  Wu  [1982].  His  work  suggests  that  instead  of  (38), 


oi  ri*  r 

Jjp:  J  d<t>  I  d6  Ip,!1  sin  (6) 


sin  (0e/2)  ■■ 


dO  |pj2  sin*  (6/2)  sin  (0) 


0C  S  30° 


Comparing  (44)  and  (38),  Sato  [1982a]  indicates  that  all 
energy  scattered  at  angles  greater  than  0e  should  be  con¬ 
sidered  lost  from  the  primary  wave  (the  first  integral  in  the 
brackets)  and  a  fraction  of  the  more  forward  scattered  energy 
should  also  be  excluded  (the  second  integral).  Wu  [1982]  ex¬ 
cluded  only  the  backscattered  energy.  To  evaluate  (44),  I  shall 
use  (26);  while  (26)  is  not  valid  near  6  »  0,  the  presence  of  the 
weighting  factor  [sin*  (0/2)  sin  (0)]  in  (44)  ensures  that  the 
integral  is  small  in  this  region.  Then 


In  f  \  dr  /,. 


from  (42)  this  gives 


- 7Wx ) 

\  W  J 


The  value  in  (47)  is  larger  than  that  in  (43)  because  a  greater 
part  of  the  scattered  energy  was  considered  lost  from  the  pri¬ 
mary  wave. 

The  results  in  (43)  and/or  (47)  may  be  compared  with  an¬ 
other  parameter,  the  backscattering  turbidity 

(48) 


(It1}  (dN\ 

9'“  in  \dr),m(l 


Hx)  =  exp  [  —  x/a]  x  >  0 


has  the  autocorrelation  given  by  (27)  if  w(x)  is  any  function 
with  a  power  spectrum  that  is  constant  with  frequency.  If  w(x) 
is  chosen  to  be  “white  noise,"  considered  to  be  a  set  of  impul¬ 
ses  infinitely  closely  spaced,  it  appears  that  (52)  may  be  con¬ 
tinuous.  even  though  it  has  an  autocorrelation  which  has  a 
nonzero  first  derivative  at  zero  lag.  However,  the  first  deriva¬ 
tive  of  (52)  is 

p'(x)  »  f  w(t)h'(x  -  r)  dr 


*  J  w(t)[0(x  -  t)-(l/a)Mx  -  r)]  dr 

-  w(x)-(l/u)p(jc)  (54) 

This  is  a  discontinuous  function,  since  w(x)  is  discontinuous. 
Thus  the  function  tdx)  is  “continuous"  but  is  not  “smooth." 
This  property  of  /4x)  is  perhaps  the  reason  that  Taylor's 
[1920]  and  Chernov’s  [1960]  proofs  fail.  The  lack  of  smooth¬ 
ness  may  explain  why  the  results  obtained  here  are  equivalent 
to  geometrical  scattering  from  obstacles,  since  in  this  case 
scattering  at  high  frequencies  is  controlled  by  the  disconti¬ 
nuity  which  is  the  boundary  of  the  obstacle.  Furthermore,  this 
discontinuity  will  scatter  waves  of  any  frequency. 

As  an  additional  limitation  on  the  theory  presented  here,  I 
note  that  there  is  a  class  of  autocorrelation  functions,  the  Von 
Karman  function  [Wu,  1982;  personal  communication,  1983] 
that  have  infinite  slopes  at.  r  =  0  for  certain  values  of  their 
parameters.  The  Von  Karman  function  is  given  by 


iGK) 


This  parameter  controls  the  amplitude  level  of  coda  waves 
relative  to  the  direct  wave,  after  correction  for  geometrical 
spreading  and  attenuation  [Aki  and  Chouet,  1975].  Aki  [1980] 
suggested  that  g„  -  g;  Sato  [19826]  presented  a  more  detailed 
calculation.  According  to  the  formula*  presented  here, 

gjg  - 1/(4*)  a  io%  (50) 

if  expression  (43)  is  used  for  g, 

gjg  a  1/(56*)  at  1%  (51) 

if  (47)  is  used  for  g. 

A  final,  important  question  has  to  do  with  the  nature  of  the 
slowness  fluctuation  idx)  that  will  lead  to  an  autocorrelation 
N(r)  that  has  a  nonzero  first  derivative  at  r  -  0.  Chernov 
[1960]  assarts  that  the  first  derivative  of  N(r)  must  tend  to 
zero  for  r-*  0  if  idx)  is  continuous,  although  a  proof  is  only 
given  for  the  one-dimensional  case.  This  is  a  special  case  of  a 
result  due  to  Taylor  [1920],  who  proved  (for  the  one- 
dimensional  case)  that  the  autocorrelation  of  a  continuous 
function  has  only  even-ordered  derivatives  that  are  nonzero  at 
zero  lag.  Pekeris  [1942]  comments  further  on  Taylor's  result. 
However,  K.  Aki  (personal  communication,  1983)  has  pointed 
out  that  the  (one-dimensional)  function 

l4x)  m  f  w(T)h(X  - 1)  dt  (52) 


where  T(m)  is  the  gamma  function  and  Km  is  'the  modified 
Bessel  function.  For  m  <  1/2.  ( dN/dr),.0~*  <x>.  Obviously,  the 
theory  presented  here  cannot  apply  to  this  case,  since  we  may 
not  write  (1 1).  For  m  *  1/2.  (55)  becomes  (27),  and  for  m  >  1/2, 
(dN/dr),m0-»0. 
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